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We study one-dimensional exclusion processes in two coupled closed rings consisting of a common 
diffusive channel and two parallel active (driven) channels. Our model displays bulk-driven phase 
transition and phase coexistence in the form of a localised domain wall (DW) in one of the active 
channels in a limit where the diffusive and driven dynamics compete. By controlling a splitting 
parameter which tunes the in-coming currents into the active channels, the system can be brought 
to a delocalisation transition, when delocalised DWs are formed in both the active channels. We 
characterise the DW fluctuations numerically. 
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I. INTRODUCTION 

Totally Asymmetric Simple Exclusion Process 
(TASEP) [l[ serves as a paradigmatic example of open 
non-equilibrium systems in one dimension (Id). Its prac- 
tical realizations include quasi-ld motion of molecular 
motors along with microtubules in intra-cellular trans- 
port [2], protein synthesis Q or motion in geometrical 
confinement, e.g., nuclear pore complex of cells 14] . In 
contrast. Symmetric Exclusion Process (SEP) [5[ is a 
typical example of Id diffusion. Well-known examples 
of SEP include diffusion through artificial crystalline 
zeolites 0. In both passive (SEP) and active (TASEP 
or TASEP-like) systems, prohibition of mutual passage 
of particles or exclusion gives rise to nontrivial collective 
effects, whose details of course depend upon whether 
the dynamics in question is TASEP or SEP. Active 
systems with open boundaries generally display spatially 
nontrivial steady state density distributions. 

In this paper, we propose a closed model that con- 
sists of two overlapping rings with a common diffu- 
sive part (SEP) and two parallel active (driven) chan- 
nels (marked Ta and Tg hereafter). In order to en- 
sure competition between driven and diffusive dynam- 
ics, we consider a particular limit of the model. Our 
principal result includes identification of a model param- 
eter 9, having values between and 1 (sec below), as 
a switch, by tuning which continuously keeping every- 
thing else unchanged (i) one may de-pin pinned domain 
walls (DW) and (ii) as crosses 1/2, a localised DW 
in one of the active channels disappears and appears in 
the other. Our model should serve as a paradigmatic 
example of localisation-delocalisation transition in a Id 
closed model with coupled diffusive and driven dynamics. 
In addition to its direct theoretical relevance, it is phe- 
nomenologically motivated by the movement of molec- 
ular motors in closed compartments 3|i dynamics 
of colloidal particles in optical traps [sj and the dynam- 
ics for multiple mRNAs competing for finite resources 
(ribosomes), where the ribosomes in turn are bounded 
by a certain trajectory and a diffusion rate outside the 



mRNA (during recycling) [1, [T3|- In particular, protein 
synthesis involves two stages: transcription of genetic 
information from DNA to messenger RNA (mRNA) by 
RNA polymerase and translation from mRNA to proteins 
through ribosome translocation. In most bacteria such 
as E.coli, translation involves three main players: the 
mRNA (genetic template), the ribosome (assembly ma- 
chinery), and aminoacyl transfer RNAs (aa-tRNAs), i.e., 
transfer RNAs "charged" with the corresponding amino 
acid. The process of translation consists of ribosomes 
moving along the mRNA without backtracking. This is 
modeled by TASEP. It is well-known that ribosomes that 
move along mRNA strand are recycled in a cell. For in- 
stance, in eukaryotic cells, after each round of protein 
synthesis, the ribosomes are released from the mRNA 
and they join the common pool of ribosomes in cytoplsm, 
where they execute diffusion and may rejoin the mRNA 
to restart protein synthesis. In our model, the SEP chan- 
nel models the "common pool" of diffusive ribosomes in 
the cytoplasm of an eukaryotic cell, which in our model 
can come back to the entry point of the TASEP lanes 
due to the feedback from the SEP channel. However, al- 
though ribosome translocation along mRNA forms phys- 
ical motivation of the present work, the analogy between 
our model and the actual biological process of ribosome 
translocation along mRNA strands is not strict due to 
various limitations of our model, as we discuss below. 
In both SEP and TASEP, each lattice site has maximum 
unit occupancy and a particle can only move to the near- 
est neighbour site (in both directions for SEP or in one 
direction only for TASEP), only if that site is empty. 
Thus the dynamics obeys the exclusion principle. For 
SEP with open boundaries, the density profile is always 
linear with the slope being determined by the bound- 
ary conditions at the two ends [sj. In contrast TASEP 
with open boundaries displays three distinct phases [ll| 
characterized by their average densities (low and high) 
and a third phase marked by a maximal current (MC). 
Unlike with open boundary conditions, individual SEP 
and TASEP dynamics with closed boundaries (say, closed 
rings) exhibit only uniform density profile in the steady 
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state due to spatial translational invariance. The rest of 
the paper is organised as follows: In Sec.|IT]we discuss our 
model in details. Then in Sec. IIIII we set up our mean- 
field theory (MFT) and discuss the steady state density 
profiles by using our MFT, complemented by extensive 
Monte-Carlo simulation (MCS) studies. In Sec.|lV]we go 
beyond MFT, and discuss domain wall fluctuations and 
delocalisation transition at 6* = 1/2. Finally, in Sec. |V] 
we summarize our results. 



or Tb is a function of 9,np,D and is spatially constant. 
Notice that our model is a variant and extension of that 
in Ref. (T^]. In particular, for = 1 or 0, Tb or Ta 
is blocked and our model explicitly reduces to that of 
Ref. [ill. Evidently, for 6* > 1/2 and 6* < 1/2 the be- 
haviour of the two channels are simply interchanged. 



III. STEADY STATE DENSITY PROFILES 



II. THE MODEL 
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Figure 1: (color online) Schematic diagram of the model; 
LJ/RJ refer to the left/right junctions. Site labels run from 
m = 1 to L from LJ to RJ for Ta,b and for SEP from RJ 
to LJ, superscripts denote site numbers for the corresponding 
channel Ta, Tb and S. 

Our proposed Id model is a closed system of two over- 
lapping rings consisting of three channels of equal number 
of sites designated by m = 1, 2, ...L, as shown in Fig. ((!]). 
Dynamics of the two channels Ta and Tb are governed 
by TASEP and the particles in the third channel S exe- 
cute SEP. Thus in S particles can hop to both direction 
with rate whereas in Ta and Tb particles can only 
hop to its right neighbour if empty with rate unity set- 
ting the time scale. At the left junction of Ta and Tb 
particles can either enter from SEP channel with rate Dd 
and D{1 —9) respectively if those sites are vacant or can 
hop to the other side with rate D. If both Ta and Tb 
try to inject a particle into SEP, then one of the TASEP 
channels {Ta or Tb) is selected randomly for injecting a 
particle to the target site i.e, the first site of SEP if it is 
empty. If Np be the total number of particles then, the 
global particle density rip = Np/2>L. In Fig. ([T]), symbols 
LJ and RJ refer to the left and right junctions in the 
model. In addition, for the purpose of clarity, the site 
number for a particular channel (i.e., Ta, Tb or S) are 
given as a superscript, e.g., at LJ , the first sites of Ta and 
Tb are denoted as T\ and Tg , respectively, and the L-th 
site of S is denoted as 5*^. Similarly for the right junction 
RJ . In this model, the three bulk parameters {9,np,D) 
control different phase transitions. Thus we observe bulk 
induced phase transitions unlike the usual TASEP with 
open boundaries which display boundary induced phase 
transition. The steady state current in each of SEP, Ta 



We use mean-field theory (MFT) together with ex- 
tensive Monte-Carlo simulation (MCS) of our model to 
obtain the steady state density profiles. In the MFT, 
the system is considered as a collection of three chan- 
nels (two TASEP and one SEP) with effective entry and 
exit rates [l3|. Once these effective rates are determined 
from the condition of constancy of particle currents, one 
may use them in conjunction with the known results 
for TASEP and SEP with open boundaries to obtain 
the density profiles here. Since an isolated TASEP in 
steady state can be in three different states, the low den- 
sity (LD), high density (HD) and maximal current (MC) 
phases, and we have two active (TASEP) channels, there 
are a number of possibilities for the overall density pro- 
file of the two active channels. In order to ensure that 
the diffusive current does not vanish in the thermody- 
namic limit (TL, see Ref. [11], see below also) we let 
diffusivity D scales with system size L and define a pa- 
rameter d ~ D/L which is the same for any arbitrary 
system size. Thus steady states of the model are to be 
parametrised by (d. Hp, 9). Let us now set the notations: 
for discrete lattice, density at a particular site m is de- 
fined as = (n™), where i = A and B refer to Ta 
and Tb, and i = S for the mean density in the SEP 
channel S. Further in MFT considering continuum limit 
the density is defined as Pi{x), where x = m/L, and in 
TL, L ^ 1, X lies in the range < a; < 1. In all our 
MFT analysis we use the continuum labelling x for the 
lattice in one dimension. Our main results are summa- 
rized in the phase diagrams (parametrised by Up and d) 
for pa{x) and pb{x), as shown in Fig. for a repre- 
sentative value oi 9 = 0.8, which are obtained by exten- 
sive Monte-Carlo simulations of our Id lattice-gas model 
and the corresponding Id MFT; see below for details. 
The Monte-Carlo simulations were realized by random 
sequential update. For any value of 9 (except 9 = 0,1, 
when one of the active channels is closed), the phase di- 
agrams of both Ta and Tb may display a combination 
of the usual LD, HD and MC phases and a region of 
co-cxistence of LD and HD phases, i.e., when the density 
profiles show localised domain walls (DW). The phase co- 
existence regions arc non- overlapping for Ta and Tb, i.e., 
they do not appear for the same values of d and Up for a 
given 9^1/2. kt 9 = 1/2, for which Ta and Tb are sta- 
tistically symmetric, the density profiles pa{x) and pb{x) 
of Ta and Tb, respectively, are naturally identical. The 
crucial difference with 9 ^ 1/2 is that the co-existence 
region now corresponds to delocalised DWs, i.e., the frac- 
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Figure 2: (color online) Phase diagram with 9 = 0.8 for densi- 
ties Pa and pB for channels Ta and Tb respectively obtained 
from our MFT (dashed and dotted lines) and MCS (circles, 
squares and triangles) analysis. Staying up to d = 1, Ta has 
four phases (all marked by sufhx A in the phase diagram), 
LD,LD-HD,HD and MC phase, where as, Tb has only three 
phases (marked by sufSx B) without any MC phase. The dot- 
ted blue and the dashed black lines correspond to the phase 
boundaries obtained for Ta, whereas the dashed-dotted or- 
ange line for Tb from MFT calculations. The red circles and 
grey triangles correspond to the phase boundaries obtained 
from MCS for Ta and the green squares for that of Tb ■ 



tion of the system size L visited by a fluctuating DW does 
not vanish in TL L — >■ cx). Such dclocalised DWs appear 
in both Ta and Tb only for = 1/2. Thus as e ^ 1/2, the 
model undergoes a delocalisation transition. The phase 
diagram for 6 = 1/2 are given in Fig. 




Figure 3: (color online) Phase diagram for both Ta and Tb 
for 6 = 1/2. There are four distinct phases; Ta and Tb, 
being statistically identical, have identical phase behaviour 
(see text). The LD-HD phases in both Ta and Tb now have 
delocalised DWs. Numerical data are shown by circles and 
triangles and the phase boundaries obtained analytically are 
shown by solid lines. 

To begin with, we denote densities at the junction site 
by a^.B = Pa,b(0) and 1 - /3^^ = Pa,b{^) for the ac- 
tive channels, and 7 = ^5(0) and 5 = Ps'(l) in the pas- 
sive SEP channel. The SEP current then takes the well- 



Js^h- 5)d. 



(1) 



Evidently, Jg remains finite in TL provided D scales with 
L linearly, else, for a fixed D the SEP current Js vanishes 
for L — >■ 00. This provides a posteriori justification for 
the scale-dependent D that we have mentioned before. 
Noting that the current in each of Ta and Tb is given by 
Ji = ~ Pi), i = A, B (assuming no boundary layer 
at z = A,B, i.e., Ta, Tb are in their LD/coexistence 
phases) and using conservation of total current at the left 
and right junctions for the individual incoming/outgoing 
currents to/from Ta and Tb from/to the SEP channel 
we obtain 

JT{0) = 5{l-a,,)eD, J'B"iO) = 5{l~a,){l-9)D. (2) 

Next, the individual outgoing currents at the sites 1 in 
Ta and Tb are (again assuming that the channels are in 
LD or coexistence phase) 



), jr(o) 



(1 



a. 



J- (3) 



Conservation of current then yields </™g(l) = "^^"iCl)- 
As expected, this holds so long as Ta, Tb are in their LD 
or coexistence phase. 

In contrast, if Ta, Tb are in their HD or coexistence 
phases the total outgoing current from Ta and Tb to the 
SEP channel is given by 



jr (1) = (1 - - 7) + (1 - - 7). 



(4) 



Conservation of total current at the right junction then 
yields (assuming Ta,Tb to be in HD or coexistence 
phases) 



J°"*(1) = /3^(1-/?J+/3^(1_/?J. 



(5) 



Further, again assuming HD or coexistence phases 
for pa.b, and separately considering the currents from 
TaMx = 1) to S{x = b) yields = [ll. This is 
corroborated by our MCS simulations (see below). This 
immediately yields that the bulk currents are equal. This 
is possible only when the bulk currents in Ta,b are con- 
trolled by RJ, i.e., Ta or Tb are both in HD or a combina- 
tion of coexistence and HD phases. Notice that the con- 
ditions obtained for pA,BiO) and pa,b{^) by using current 
conservations at the respective sites do not hold simulta- 
neously, unless Ta or Tb are in coexistence phases, such 
that there are no boundary layers at a; = 0, 1 with pa,b 
being piecewise continuous. Having defined effective en- 
try and exit rates (valid separately for LD/coexistence 
or HD/coexistence phases in Ta,b) for the active chan- 
nels, we can now apply the known results of TASEP here. 
One obtains the low (high) density phases in the peri- 
odic system equally and are characterized by a uniform 
density below (above) 1/2 and a boundary layer at the 
right (left). However for ^ = /3^b, the boundaries 
are matched by a piecewise constant density profile with 
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an intervening DW. For TASEP with open boundaries, 
particle entry and exit events are uncorrelated, and as 
a result, the DW is delocalised and undergoes random 
walks covering the entire span of the system in the long 
time limit. However in the present model, as in Ref. [l^ . 
entry and exit of particles are not uncorrelated; they get 
correlated by the fact that the ends of the active chan- 
nels arc connected by the passive channel. Consequently, 
as revealed by our Monte Carlo simulation studies, we 
find localised DW in the active channels, which is sim- 
ilar to Ref. [Hj. However, rather surprisingly for the 
special case oi 9 = 1/2, i.e., when each of the active 
channels carry equal current on average, we obtain de- 
localised DWs in both channels. Our MFT formulated 
above may now be used to analyse the density profiles in 
the different channels of the model quantitatively. 



A. DW in one active channel and LD in other 

First, let us consider a situation when there is a DW 
in one of the active channels (say Ta with 9 > 1/2) and 
the other active channel (Is) is in the LD phase with a 
uniform density (within MFT neglecting any bound- 
ary layer). Following the phenomenology of TASEP with 
open boundaries, we set = as a requirement of a 
DW in Ta- Possibilities of simultaneous DWs in Ta and 
Tb will be discussed later. Within MFT, pA may be rep- 
resented by a Heaviside function that connects the two 
regions of constant density through a localized DW at 
position x^, (say) as. 



(6) 



Since Tb is assumed to be in the LD region, density pB (x) 
can be written as, 



PBix) = 



(7) 



neglecting the boundary layer at the right boundary. For 
the SEP channel the linear density distribution gives. 



p{x) = (5 -I- (7 - (5)a 



(8) 



Further, the particle number conservation can be ex- 
pressed as 

3np = / dx[ p{x) + pa{x) + pb{x) ], (9) 
Jo 

following the conditions as above and disregarding the 
discontinuities at the right boundaries. Again from 
Eq.(l2]) and Eq.^ we have. 



1 



- 1 = a^q, 



(10) 



where q = (1/6*- 1). Now by solving Eq. ^ in TL {S 0) 
we get 



1 + § - Srip - a^(l - q) 
l-2a. ' 



(11) 



Again Eqs. (0]), (O, pO)) and the relations = and 
= 1 — yield for 7 as. 
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alil + q^)~2a^+l 
l-a,{l-q) 



(12) 



In TL (5 — and 7 = Jjd. Since the model considered 
here is closed, J, = J^"*(0)-h J]^"*(0). Again from Eq. © 
we have, 



"A(l-"A) + (l-g"A)g«^ 
7 



(13) 



Hence, the position of the DW depends on the two con- 
trol parameters rip and d for a given 9. When the DW in 
Ta is localised within the system (0 < a;^ < 1), it con- 
nects the LD and HD phases of Ta through a phase of 
coexistence (LD-HD). The boundaries between the LD, 
LD-HD phases and LD-HD, HD phases of Ta are ob- 
tained by setting x^, = and x^ = 1 respectively. Set- 
ting x^ = from Eq. pT|) we get a quadratic equation in 
q;^ . For a particular value of Up the feasible values will 
be (0 < < 1/2). Now putting that in Eq. (HH) we 
get the corresponding d value. Thus we get the bound- 
ary between LD and LD-HD coexistence phase in the 
{rip, (i)-plane. Similar exercise for x^ = 1 gives the right 
boundary between LD-HD and HD phase. See Fig ([2]) 
for details. In Fig ([2|) the phases of channel Ta and Tb 
are spanned by c? < 1 and Up < 1. The phase diagrams 
for Ta and Tb are drawn corresponding to a situation 
when Ta displays a variety of phases (LD, LD-HD, HD 
AND MC), while Tb remains in its LD phase. For d > 1, 
this part of the phase diagram remains qualitatively un- 
changed, with the phase marked as MCa, LDb (i.e., Ta 
in MC and Tb in LD) should expand to a larger area. 
Similarly, the phases of Tb are shown when channel Ta 
remains in its HD phase. The condition for the latter is 
mathematically given by pa{x) = 1 — d/2 in the bulk. 
For 1 < 0? < 2,pa{x) < 1/2 and hence Ta is no longer 
in its HD phase. Hence, we do not consider the d > 1 
region while presenting our phase diagram. A DW in 
Ta obtained from our MCS studies are shown in Fig. (jH 
top) with 9 = 0.80, Up = 0.40 and d = 0.15. We have 
taken L ~ 100 and 200 for determining the position of 
the domain walls and phase diagram. We do not find any 
significant dependence of our results on L. 



B. DW in one active channel and HD in other 

Let us now consider the case when there is a DW in 
Tb, and Ta is in HD phase for 9 > 1/2, thus having a 
boundary wall at the left end. As discussed above, within 
MFT, pb may be represented by Heaviside 6*- function as 



PBix) 



+ eix~x^){i 



(14) 



and Ta is in HD phase having a uniform density of (1 — 
f5^), neglecting the boundary layer. To have a DW in Tb 
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Figure 4: (color online) (top) For L = 200, 9 = 0.8, Up = 0.40 
and d = 0.15, Tb is in LD phase and DW appears in Ta at 

0.19. (bottom) 
Uj, = u.oo ana a = 0.20, Ta is in HD 



x-^ - 0.21, whereas from MFT Eq.(ITTl) 
For L = 200, 9 = 0.8, n„ = 0.65 and d 



phase and DW appears in Tb at ~ 0.38, whereas from 



the domain walls obtained from MFT. 



wc must have 
get 

7 ^ 



/3^. Now from Eqs. (HJ and (O, we 



(1-/3J^ + (1-/3J^ 



(15) 



Hence Eq. and /S^ yield 7 = (1 — a^) . In 

TL 7 = Js/d, and hence, = (i/2. Again by using the 
particle number conservation and as Pa — Pb — '^b ^'^'^ 
Pa{x) = 1 — (i/2 wc obtain, 



= 1 - 



3n„ — 



(16) 



From the above expression we get the boundaries between 
the LD, LD-HD {x^ = 0) and LD-HD, HD phases {x^ = 
1) of Tb. a DW in Tb, obtained in our MCS studies, is 
given in Fig. bottom). There is a crucial difference 
between the DWs in Ta and Tb: The LD part of the DW 
in Ta has density , different from the density of Tb 
(fully in LD), thus pa{x) has no overlap with pb{x)- In 
contrast, pa{x) (fully in HD) overlaps with pb{x) in the 
HD part of the DW in Tb- This is due to = (3^ and 
is clearly visible in Fig. 



C. Delocalised domain wall at 9 = 1/2 

Let us now carefully consider the properties for 6 = 
1/2, when both the active channels are symmetric and 
statistically identical. Thus, if = then automati- 
cally = Pb- Hence, if Ta has a DW, Tb too will have 
a DW, or is in its LD-HD (co-existence) phase as well, 
such that its density may be represented by a Heaviside 
^-function that connects the two regions of constant den- 
sity through a localized DW at x^ (say). Hence, pa{x) 
and pb{x) are given by the expressions ^ and (|14p re- 
spectively. As both Ta and Tb show DWs, thus = Pa 
and = Pb^ ^^'^ ^o^' ^ = from Eq.(I21) we have, 
= P^ =: = Pb = ct (say). Again, from Eq.(I2]) and 
Eq.dU we have 7 = 1 — a. Therefore, in TL, J — > and 
thus a = d/2. Now by particle number conservation we 
have. 



3np = 1/2 + 3d/4 + (2 - 4 - x^){l - d). (17) 



Thus for 6* = 1/2 we get a relation given by Eq. fTTl) 
between x^ and x^ for a particular value of Up and d. 
In other words x^ and x^ are not uniquely determined. 
In this case both DWs are delocalizcd and perform ran- 
dom walk along the active channels. As Ta and Tb 
are identical, the condition for DW is satisfied for both 
the channels simultaneously. Thus, x^ = x^ = and 
= x^ = I in Eq. (fT7|) can give the boundaries of LD- 
HD phase in both channels with the LD and HD phases 
respectively. In Fig. ([3]) we have shown this mean-field 
result as well as that obtained from MCS which shows 
distinct four phases with phase boundaries for both chan- 
nels arc identical. However, symmetry between the two 
active channels dictate that the long-time averaged po- 




Figure 5: (color online) Plots of p = pA, pB versus site (x) for 
9=1/2 with Up = 0.52 , d = 0.19 and L = 2 x 10^ Clearly, 
both PA and ps display overlapping delocalised DWs. 



sitions of the DWs (equivalently, the long-time average 
density profiles pa{x) and pb{x)) in Ta and Tb arc iden- 
tical, a fact verified by our MCS, is displayed in Fig. ([S]). 
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D. Both the active channels in MC phase 

Lastly, we consider the possibility of the MC phases in 
the active channels. Let us first consider the conditions 
for obtaining MC phases in both the active channels. 
Condition for MC in isolated TASEPs are = ps = 1/2 
and Ja = Jb = 1/4, and this happens when all the 
boundary densities aA,C(B, I3a, I^b > 1/2. Furthermore, 
Ta, Tb have boundary layers at both the ends. This 
precludes usage of Eq. ([3]) to determine the boundary 
densities. Using = 1/4 = together with Eq. ([2]) 
(assuming no density discontinuities between ps(l) and 
PA,B{0)),weRndS{l-aA)eD = 1/4, 6{l-aB){l~0)D = 
1/4. Since aA,B > 1/2 



5 > 



1 



1 



2dLe' 2dL{l - I 



(18) 



Using similar considerations at RJ, and again assum- 
ing no density discontinuity between ps(0) and jOA.s(l), 
which means 1 — 7 = (3a = Pb, together with > 
1/2, Eq.[T]and Js = 1/2 we have, 



5 < 



1 

2d' 



From particle conservation we have, 



3n, = ^+l 



1 



(19) 



(20) 



Eqs. (UHl), (HH) and (gOl) yield boundaries of the MC phase 
with LD and HD phase as, d = l/{12np — 4) and d ~ 
1/(6 — 12rip) which indicates the presence of such phase 
for d > 1 and bounded by the two above mentioned lines. 
Thus, the demarcating lines are independent of 6. They 
are shown in Fig. ([3]). We now consider the case when one 
of the channels (say, Ta) is in MC phase and the other 
one (Tb) in the LD phase. Therefore, we have Ja = 1/4, 
PA 1/2 and a^,/3^ > 1/2. From Eqn. (0) and © we 
have a 4 = Dd9. Again using the MC phase condition 
q;^,/3a > 1/2 we have (arguing as before). 



X 1 X 1 

'^>m^''^<2-T- 



(21) 



The maximal current condition gives Jg = (1/4 + ^7/2 
q^/4), then from particle conservation we have. 



Srir, 



S + 



2d 



(22) 



The two inequalities ((2T)) together with Ea. (l22t then 
yield boundaries of the MC phase with the LD and 
HD phase respectively as, d = Js/{6np — q — 1) and 
d = Js/{q + 2 — 6np). In Fig. ^ we have shown the 
MC phase boundaries. Our MCS studies also reveal a 
small MC phase within the region obtained from MET. 
Not surprisingly, for 6* = 1 and 9 = 0, the MC phase 
regions obtained from our MET match exactly with that 



of Ref. [1^. In addition, one may argue that the co- 
existence of Ta in HD and Tg in MC is not possible. 
For Ta to be in HD phase 1 - > 1/2 or < 1/2. 
Again = < 1/2, so long as Ta.b are in HD or 
coexistence phases. But condition for MC phase in Tb 
is a^jP^ > 1/2. Thus, an MC phase in Tb (when Ta in 
HD) is not allowed. 



IV. DOMAIN WALL FLUCTUATIONS AND 
DELOCALISATION TRANSITION 



Until now we have considered the MET for the model, 
where all fluctuations are neglected. However, the DWs 
fluctuate about their ME DW (mean) positions or 
x^. We have studied these fluctuations numerically and 
characterise them by measuring the scaling of the fluc- 
tuations with 6 and system size L. In particular as 
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Figure 6: (color online) Log-log plot of DW width a versus 
16* — 6'c| for L = 200, Up = 0.40 and d = 0.15 with the exponent 
/3 = 3/4 (see text). 
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Figure 7: (color online) Log-log plot of a versus L for Up = 
0.40 and d = 0.15 at an off-critical point: 9 ^ 9c- The line 
has a slope 1/2; thus, a ~ l/^/T, (see text). 

9 ^ 9c = 1/2 from above or below, the localised DW 
in Ta or Tb shows a delocalisation transition at which 
DW fluctuations diverge. The width <t of the distribu- 
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tion of DW fluctuations can be obtained by fitting the 
density profile in the vicinity of the domain wall by the 
function (P • erf[{x — Q)/a] + R) [HI, with the param- 
eters P, Q, R, a. We find a to diverge with a power law 
dependence on {9 — 6c) as, 

a^{ec-0)-P, (23) 

with /3 = 3/4 obtained from our MCS studies as shown 
in Fig. ([5]). In contrast, at an off-critical point, i.e., for 
9 7^ 1/2, DW fluctuations are finite and vanish in TL 
L — )• oo as L~^/^; wc have shown this in Fig. For 
investigating the variation of domain wall width with L 
for various values of we have taken L — 32, 64, 128, 256 
and 512. 



V. SUMMARY AND OUTLOOK 

Analytical and numerical studies of our model am- 
ply illustrate the underlying rich phase behaviour, in- 
cluding a delocalisation transition, unexpected in a sys- 
tem without boundaries. While boundary-induced phase 
transitions including delocalisation transitions have been 
observed in several open systems with exclusion pro- 
cesses together with spatially nontrivial steady state den- 
sities [l5l - [l7t . analogous studies on bulk closed systems 
are less studied so far. The competition between the dif- 
fusive and driven dynamics, and the division of the SEP 
current into two parallel TASEP currents are crucial to 
the macroscopic behaviour we obtained. The latter is 
controlled by a parameter 9, which is a tuning parameter 
in the model. The most striking feature in our work vis- 
a-vis the results in Ref. [l2j is the possible existence of 
a delocalisation transition and correspondingly the for- 
mation of DWs in both Ta and simultaneously at a 
special value 9 = 1/2. In contrast to the DWs formed 
either in Ta or Tb (but not simultaneously in both) for 
9 ^ 1/2, as found in Ref. 12] as well as in the present 
work, the DWs at 9 = 1/2 are no longer pinned to a 
fixed point in the lattice with vanishing fluctuations in 
the thermodynamic limit. Instead they delocalise and 
have position fluctuations that do not vanish in the ther- 
modynamic limit. Thus the parameter 9 in our model 
appears as a tuning parameter or a switch, which can 
be used to control the nature of domain wall fluctuations 
(localised/delocalised). In addition for 9 ^ 1/2, the value 
of 9 can be tuned to make the DW appear or disappear 
in one of the active channels. There is no analogue of 
these in the study of Ref. [l^] ■ While we have considered 
only two TASEP channels, many more may be added 
and studied systematically as above. Recalling protein 
synthesis by ribosomes along mRNA strands as one of 
the phenomenological motivation for our model, it may 



be noted that several mRNAs compete for same resources 
(ribosomes) in a cell. Thus a systematic study of multiple 
TASEP channels connected in parallel with a single SEP 
channel would be useful. The failure of the traditional 
MET calls for further analysis by means of more sophis- 
ticated analytical techniques, e.g., Bethe ansatz [131 oi' 
density matrix renormalisation group , which are be- 
yond the scope of the present work. From the point of 
view of noncquilibrium statistical mechanics, our model 
belongs to the class of models lacking translation invari- 
ance and without boundaries that displays a phase tran- 
sition (in the form of a delocalisation transition). Our 
model may be extended in several new directions, e.g., 
again motivating by ribosome movements along mRNA, 
one may in our model consider particle exchanges be- 
tween Ta and Tb, or between one of the active channels 
and passive channel (representing ribosome attachments 
or detachments), allow defects along the active channels 
(representing defects in the mRNA), introduce a second 
control parameter at the exit ends of Ta and Tb that con- 
trols the relative outgoing currents to SEP and unequal 
hopping rates in Ta and Tb- These will be considered 
elsewhere. We close this work with a note of caution: As 
mentioned in the beginning, despite some similarities our 
model cannot be directly used for quantitative descrip- 
tions of ribosome translocations along mRNA strands 
due to its limiations. First of all, ribosome diffusion takes 
place inside a cell, which, alothough geometrically con- 
fined, has a three-dimensional (3d) structure, as opposed 
to our Id diffusive model for it. Secondly, the description 
of a ribosome as a single unit (i.e., a point particle here) 
is also questionable, for it gets released from an mRNA 
by falling apart into different subunits, a feature not pos- 
sible to capture in our simplified description here. Nev- 
ertheless, our work provides some clues about the actual 
biological system, e.g., the crucial role of particle num- 
ber conservation in determining the nature of the steady 
states. We expect that more realistic theoretical descrip- 
tions of ribosome translocation and detachment should 
have some of the basic features of our model in-built into 
it. 
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